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DEFINITIONS  AND  STANDARD  CONDITIONS 


Symbols  used  in  this  report  are  tabulated  below.  Values  marked  with  an  asterisk 
(^')  refer  to  the  values  of  parameters  under  standard  conditions.  These  quantities  are 
altered  in  one  of  the  graphs.  Values  marked  with  a  plus  sign  (+)  refer  to  the  values  of 
parameters  altered  for  the  purpose  of  comparison  with  experiment. 


A  transverse  cross-sectional  area  of  graphite  lattice  available  for  heat 
flow  =  378-  cm^ 

Cg^  specific  heat  of  air  in  channel  =  0.24  cal/g-<^C* 

Cj  power  transfer  coefficient  by  conduction  from  metal  to  carbon  = 

0.0095  cal/cm-sec -OC* 

power  transfer  coefficient  by  radiation  from  metal  to  carbon  = 

1.1  X  10“^^  cal/cm-sec -OA"^* 

D  4  X  hydraulic  radius  of  air  channels  =  2.27  cm 
e  distance  of  flux  maximum  from  gap  face  of  pile  =  -3.5  cm*+ 

f  distance  of  plenum  end  of  fuel  rod  from  flux  maximum  =  339  cm*+ 

h  heat  transfer  coefficient  from  metal  or  carbon  to  air  in  channel  = 

0.0038  cal/cm^- sec-^C*"^ 

h^  heat  transfer  coefficient  from  carbon  to  air  in  gap  =  0.00093 
cal/cm^  -  sec 

K  cross-sectional  area  available  for  air  flow  in  channel  =  27.9  cm^ 
k  thermal  conductivity  of  graphite  ^  =  0.17  cal/cm -sec -OC*+ 
kg^  thermal  conductivity  of  air  =  7.6  x  10”^  cal/cm-sec-oC 
L  axial  distance  from  flux  maximum  to  point  where  flux  extrapolates  to 
zero  =  381  cm*"^ 

Li  distance  from  center  of  pile  to  edge  of  orificing 

Li'  distance  fromi  center  of  pile  to  edge  of  loaded  lattices,  i^^  direction 

N  number  of  channels  loaded  =  1369“^ 

P  power  produced  by  reactor  =  28  megawatts'*" 
perimeter  of  the  carbon  =  21.3  cm 

P^  perimeter  of  the  metallic  fuel  elements  including  the  fins  =  27.8  cm 
Q  power  generated  in  central  half-channel  =  6040  cal/sec"*" 

Q'  power  generated  per  unit  length  in  central  half-channel 
R,S  coefficients  of  the  two  particular  integrals  of  a  certain  linear,  inhomo- 
genous,  first-order  differential  equation 
s  transverse  extrapolation  distance  (=  61  cm  for  comparison  with  experi¬ 
ment  only) 

T^  temperature  of  air  in  channels 
T^  temperature  of  the  carbon 
Tm  temperature  of  the  metallic  fuel  element 
U  See  equation  (24)  for  its  definition 
V  See  equation  (22)  for  its  definition 

W  mass  rate  of  air  flow  thru  central  channel  =138  g/sec*+ 

X  horizontal  positional  coordinate  in  a  direction  transverse  to  that  of  the 
fuel  elements 

y  vertical  positional  coordinate  in  a  direction  transverse  to  that  of  the 
fuel  elements 


^R.  C. 


Garth  and  V,  C, 


Sailor, 


“Thermal  Conductivity  of  Graphite,  “  BNL-69, 


November  28,  1949, 
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axial  positional  coordinate  in  direction  parallel  to  that  of  the  fuel  element 
fraction  of  all  heat  generated  in  the  carbon  =  0.05*+ 
viscosity  of  air  =  225  micropoises 

two  roots  of  a  quadratic  equation.  These  roots  are  the  coefficients  of  z 
in  the  exponentials  corresponding  to  the  solutions  of  a  certain  linear, 
homogenous,  first-order  differential  equation 
tt/ZL  =  4.12  X  10-Vcm 


PREFACE 


The  axial  temperature  distribution  is  calculated  in  a  pile  with  a  central  air  gap. 
The  effect  of  various  geometrical  and  thermal  constants  on  this  distribution  is  pre¬ 
sented  in  Figures  4  through  10,  and  a  tabular  summary  is  presented  in  Figure  11. 

The  coefficient  of  heat  transfer  to  the  air  and  the  specific  heat  of  the  air  are  the  most 
important  parameters.  The  fraction  of  heat  generated  in,  and  the  heat  conducted  or 
radiated  to,  the  graphite  and  the  consideration  of  the  actual  flux  distribution  are  of 
secondary  importance.  The  thermal  conductivity  of  graphite  and  the  heat  transferred 
to  the  gap  are  of  negligible  importance. 

A  simple  theory  is  compared  with  a  more  elaborate  theory  in  Figure  3.  The  air 
temperature  in  all  cases  can  be  computed  just  as  well  from  the  simple  theory  as  from 
the  more  elaborate  one.  The  temperatures  of  the  fuel  elements,  the  air  and  the 
graphite  at  the  plenum  end  are  all  very  nearly  equal,  and  the  elaborate  and  simple 
theories  lead  to  essentially  the  same  result  for  this  temperature.  For  the  Brookhaven 
reactor  at  28-Mw  power,  the  elaborate  theory  predicts  the  fuel  element  temperature 
at  the  gap  to  be  some  ZO^C  lower  and  the  graphite  temperature  to  be  some  40®C  higher 
than  does  the  simple  theory.  These  differences  gradually  vanish  as  one  moves  along 
toward  the  plenum  chamber. 

Finally,  in  Figure  12,  an  attempt  is  made  to  compare  theory  and  experiment.  The 
agreement  is  poor  because  of  experimental  difficulties  and  a  theoretical  difficulty  as 
explained  in  the  final  section.  Numerical  Results  and  Conclusions. 
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AXIAL  HEAT  DISTRIBUTION  IN  BROOKHAVEN  REACTOR 

EFFECT  OF  VARIOUS  GEOMETRICAL  AND  THERMAL  CONSTANTS  AND 
COMPARISON  WITH  EXPERIMENT 


Introduction  and  Approximations 


The  primary  purpose  of  this  report  is  to  calculate  the  temperature  distribution 
along  the  air  cooling  channels  in  the  Brookhaven  reactor.  Not  only  is  the  temperature 
distribution  needed  in  order  to  design  a  reactor,  but  it  is  requisite  for  the  evaluation 
of  the  temperature  coefficient,  xenon  poisoning,  stability  under  flashes,  etc. 

The  secondary  purpose  of  this  study  is  to  determine  the  influence  on  this  dis¬ 
tribution  of  various  geometrical  and  thermal  constants  so  that  simplifications  can  be 
made  in  future  calculations  by  neglecting  unimportant  considerations. 

A  third  purpose  is  to  compare  a  simple  theory  and  a  more  elaborate  theory. 
Fourthly,  the  elaborate  theory  is  to  be  compared  with  experiment.  Unfortunately,  a 
valid  comparison  is  almost  impossible  because  of  experimental  reasons  and  because 
the  theory  in  at  least  one  respect  is  inadequate.  The  reasons  are  presented  in  the 
section,  Numerical  Results  and  Conclusions. 

The  Brookhaven  reactor  consists  of  a  graphite  cube  split  down  a  central  plane. 

Air  flows  into  this  split  or  gap  and  out  of  the  air  channels,  as  shown  in  Figure  1.  The 
air  channels  of  the  Brookhaven  reactor  have  their  flow  restricted  or  orificed  in  such  a 
a  way  that  the  mass  rate  of  flow  through  any  channel  is  proportional  to  the  heat  energy 
produced  in  that  channel.  This  fact  introduces  the  simplification  that  the  reactor  is, 2 
on  the  average,  approximately  isothermal  in  any  plane  parallel  to  the  gap.  We  may, 
therefore,  restrict  ourselves  to  a  consideration  of  the  central  channel  and  its  associa¬ 
ted  lattice.  Were  it  not  for  the  orificing,  this  channel  would  be  the  hottest.  We  see, 
also,  that  heat  flow  in  directions  parallel  to  the  gap  may  be  neglected  in  a  first  ap¬ 
proximation.  We  consider,  therefore,  only  the  flow  of  heat  axially  along  the  channels 
and  the  interchange  of  heat  into  and  out  of  the  ends  of  the  graphite  lattice. 

The  problem  of  the  temperature  distribution  within  the  fuel  element  itself  has 
been  considered  elsewhere  and  will  not  be  treated  here.^  The  temperature  of  the  fuel 
element  will  be  characterized  by  its  surface  temperature  at  any  axial  point. 

The  radiation  of  heat  from  the  fuel  elements  is  proportional  to  the  fourth  power 
of  the  temperature.  Because  this  nonlinear  dependence  introduces  considerable 
mathematical  complications,  we  linearize  this  dependence  by  considering  just  the 
first  two  terms  in  a  Taylor  series  expansion.  It  is  to  be  noted  that  only  about  20%  of 
all  the  heat  is  transferred  by  radiation  (to  the  graphite),  most  of  the  remainder  being 


^It  is  to  be  noted  that  the  heat  transfer  coefficient  is  proportional  to  the  mass  rate 
of  air  flow  to  the  0.8  power.  Hence,  the  periphery  is  slightly  cooler. 

^M.  Clark,  Jr.,  “Temperature  of  and  Stresses  in  Cylindrical  Fuel  Elements  During  Pile 
Flashes,  “  BNL-86. 
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transferred  directly  to  the  air.  The  linearization  approximation  will  not  introduce  a 
serious  error. 

Heat  may  be  transferred  to  the  graphite  from  the  fuel  elements  by  radiation  and 
by  conduction,  and  from  the  air  by  conduction.  Further,  heat  may  be  generated  in  the 
graphite  by  the  slowing  down  of  neutrons  and  by  the  absorption  of  gamma  rays. 

The  air  in  the  channels  is  so  highly  turbulent  that  it  may  be  Considered  thoroughly 
mixed  across  any  cross  section  transverse  to  the  channel.  The  air  may  at  any  axial 
point  be  characterized  by  a  single  temperature. 

We  shall  assume  the  power  generated  varies  co sinusoidally  from  the  center  of 
the  gap."^ 


Simple  Method  of  Calculation 


Before  formulating  the  more  accurate  method  of  calculation,  let  us  consider  a 
simple  first  approximation.  We  neglect  all  forms  of  heat  transfer  except  the  conduc¬ 
tion  from  the  fuel  elements  to  the  air  stream  and  from  the  carbon  to  the  air  stream. 
Three  statements  serve  to  formulate  the  equations: 

1.  The  power  generated  by  the  metal  is  equal  to  that  carried  away  by  the  air: 


TTQ  7TZ 

(1-a) - cos -  =  hP^  (T 

2L  2L  m  ^  m 


(1) 


where 


a  =  fraction  of  all  energy  of  fission  released  in  carbon 
Q  =  power  generated  in  central  half-channel  at  flux  maximum 
L  =  half-length  of  extrapolated  pile  along  fuel  element  direction 
z  =  distance  of  point  of  interest  along  a  fuel  element  measured  from 
the  center  of  the  reactor 

h  =  power  transfer  coefficient  to  the  air  stream  from  either  metal  or 
carbon 

Pm  =  perimeter  of  metal  fuel  element  including  fins 
Tm  temperature  of  the  metal 
T^  =  temperature  of  the  air 

2.  The  power  generated  in  the  carbon  is  equal  to  the  power  carried  away  by 

air: 


77Q  77Z 

a -  COS  -  =  hP  (T  “T  ) , 

2L  2L  c  V  c  a/ 


where 


^J,  Chernick,  I.  Kaplan,  J.  Kunstadter,  V,  Sailor  and  C*  Williams,  “An  Experimental 
and  Theoretical  Study  of  the  Subcritical  BNL  Reactor,**  BNL-60* 


(2) 
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=  perimeter  of  carbon 
=  temperature  of  carbon 


3.  The  power  absorbed  by  the  air  equals  the  power  absorbed  from  the  car¬ 
bon  and  from  the  metal: 


CaW  —  =  (Tc-Ta)  +  (T^-T^), 


(3) 


whe  re 


=  specific  heat  of  the  air 
W  =  mass  flow  rate  of  air  in  the  central  channel 

The  above  equations  are  readily  solved  by  finding  one  temperature  in  terms  of 
the  others  and  subsequently  eliminating  that  temperature  in  the  remaining  equations. 
The  results  are: 


0  7TZ  (4) 

=  -  sin - 

®  C  W  2L 

a 


T 


m 


(l~a)77Q 

2hP^L 

m 


7TZ 

COS  -  + 

2L 


(5) 


T 


c 


arrQ 

2hP^L 


77Z 

COS -  +  T, 

2L  ' 


(6) 


The  zero  of  the  temperature  scale  is  taken  to  be  the  inlet  temperature  of  the  air. 


The  mass  flow  rate  of  air  in  the  central  channel  is  easy  to  determine.  Because 
of  the  orificing,  the  air  flow  in  a  channel  varies  as  cos  (ttx^/ZLj)  cos  ('rTy'/2L2)»  where 
the  channel  is  located  at  (x^,  y^).  Hence,  the  total  flow  of  air  is: 


2W  S  cos 

i  .  j 


(7) 


The  summation  may  be  replaced  by  an  integration;  the  total  air  flow  is: 


2WN 

4L,L2 


(8) 


8NW 


2 

77 


(9) 


where  N  =  number  of  channels  loaded  into  a  rectangular  array,  2Lj^  by  in 

size.  Since  the  total  mass  flow  rate  of  air  is  known,  it  is  simple  to  determine  W 
from  equation  (9). 
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The  power  transfer  coefficient  is  determined 


5 


from  the  relation; 


hD 

/  DW 

-  =  0.023 

1 

\  J^K 

(10) 


where 


=  thermal  conductivity  of  air 
D  =  4  X  hydraulic  radius 
jjL  =  viscosity  of  air 

K  =  cross-sectional  area  available  for  air  flow 
Cg^  =  specific  heat  of  air  at  constant  pressure 


Finally,  Q  is  determined  from  the  total  power  P; 


P 


2  2 
i  .  j 


cos 


It 


8QN 


(11) 


(cf.  equation  (9)). 


Elaborate  Method  of  Calculating  Temperature  Distribution 


We  turn  now  to  a  more  elaborate  method  of  calculation.  Three  statements  of 
energy  conservation  suffice  to  establish  the  equations  of  heat  flow: 

1.  The  power  made  in  the  fuel  element  equals  the  power  carried  away  by 
conduction  to  the  carbon  plus  that  carried  away  by  radiation  to  the  carbon  plus  that 
carried  away  by  convection  to  the  air: 


/  TTZ 

(1-a)  Q  cos 

2L 


Cl  (T„-Tc)  +  C, 


+  hP 


(T„-T„), 


(12) 


where 

Q'  =  power  generated  at  flux  maximum  per  unit  length  of  fuel  element 
=  irQ/2L 

z  =  axial  distance  to  point  of  interest  from  maximum  of  neutron  flux 
^W.H.  McAdams,  Heat  Transmission.  2nd  ed. ,  McGraw-Hill  Book  Co.,  Inc.,  1942,  p.  168. 
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=  power  conduction  coefficient  from  metal  to  carbon 
C2  =  power  radiation  coefficient  from  metal  to  carbon 

2.  The  power  produced  in  an  axial  section  of  the  carbon  equals  the  change 
in  power  transported  parallel  to  a  channel  plus  the  power  conducted  to  the  metal  plus 
the  power  lost  by  radiation  to  the  metal  plus  the  power  carried  away  by  the  air: 


aQ  cos -  —  kA  - r 

2L 


-  C,  (T„-T^)  -  Cj  (T„-T^)  +  hP^  (T^-T3) 


(13) 


where 


k  =  conductivity  of  carbon 

A  =  transverse  cross-sectional  area  of  carbon  in  one  lattice 

'3.  The  increase  in  temperature  of  the  air  multiplied  by  the  heat  capacity 
and  the  mass  rate  of  air  flow  equals  the  power  given  to  the  air  from  the  metal  plus 
the  power  given  to  it  from  the  carbon: 

dT 

C^w-^  =  '’Pm  ^Tm-Ta)  + 


The  boundary  conditions  are: 

1.  The‘ power  flowing  axially  from  the  graphite  to  the  air  in  the  gap  is  pro¬ 
portional  to  the  temperature  difference  between  the  carbon  and  the  air  at  that  point. 
At  z  =  -|  e  I 


k 


h' 


(15) 


where 


|e|=  the  distance  from  the  maximum  of  the  neutron  flux  to  the  gap  face 
(see  Figure  2); 

h'=  the  heat  transfer  coefficient  to  the  air  in  the  gap  from  the  graphite 

2.  The  power  flowing  axially  from  the  graphite  to  the  air  in  the  plenum 
chamber  is  proportional  to  the  temperature  difference  between  the  carbon  and  the 
air  at  that  point.  However,  the  temperature  difference  is  small  and  the  power  transfer 
coefficient  is  small,  because  of  the  low  transverse  9.ir  velocity.  Hence,  we  may  neg¬ 
lect  the  power  transfer  to  the  plenum  from  the  face  of  the  graphite.  A  further  ap¬ 
proximation  is  in  order.  The  fuel  rods  end  at  a*  point  short  of  the  pile  face,  and  the 
flux  extrapolates  to  zero  at  a  distance  intermediate  between  the  ends  of  the  rods  and 
the  plenum  chamber  face  of  the  pile.  Since  much  of  the  heat  produced  in  the  carbon 
is  produced  by  heat  radiation  from  the  fuel  elements,  by  slowing  down  neutrons,  and 
by  stopping  gamma  rays,  we  take  the  point  where  the  flux  extrapolates  to  zero  as  the 
point  to  apply  the  zero  power  transfer  condition.  At  z  =  L, 


0. 
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(16) 


d  z 


Any  other  treatment  would  require  a  complicated  analysis  in  several  dimensions  and 
would  make  essentially  no  difference. 

The  linearized  equations  become: 


(1-a)  o'  cos 


(T  4-  T;;  (T^~T^)  +  hP^  (T  -T„) 

'in  C'  2  m'm  c-''  rn'm  b.' 


(17) 


aQ'  cos  (’’n.-Tc)  +  »’Pc  (Tc-Ta)  (18) 


dTa 


(19) 


These  equations  are  solved  by  the  tedious  process  of  systematically  solving  for 
one  temperature  in  terms  of  the  others  and  of  eliminating  that  temperature  from  all 
other  equations.  We  may  arrive  at  an  inhomogeneous,  third  order,  linear  differential 
equation  for  the  air  temperature.  This  equation  is  solved  by' the  usual  methods.  The 
complete  solution  consists  of  a  solution  of  the  homogeneous  equation  which  involves 
a  positive  and  negative  exponential,  plus  a  solution  of  the  inhomogeneous  equation 
consisting  of  a  linear  combination  of  sin  a;z  and  cos  wz  {co  =  tt/ZL).  The  two  roots  of 
the  homogeneous  equation  are: 


-kA  ±  \/(kA)  ‘  +  4C^WV 


1.2- 


2V 


where 


V 


kAC^W  (C,+4C2T^+hP^) 


hP„  (C,+4C,  T=)  +  hP^  (C,+4C,T>P„) 


(20) 

(21) 


(22) 


The  two  roots  tt^ 


and 


provide  the  coefficient  of  z  in  the  two  exponentials. 


The  coefficients  of  the  two  terms  in  sin  coz  and  cos  coz,  corresponding  to  the 
inhomogenous  solution  are  determined  by  substitiution  of  the  solution  into  the  differ¬ 
ential  equation.  The  coefficient  S  of  the  cosine  term  is: 


S 


U  +C  W) 


(Va)*+C.W)*  +  (kA)*  co^ 


(23) 


where 
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U  = 


1  + 


hP^  (l-a) 


hP.  CC,+4C3T^  +  hP^  (c,+4C2T>P,) 


(24) 


The  coefficient,  R,  of  the  sine  term  is: 


R  = 


-kAa;u 


(VW^+C.W)^  +  (kAo;)^ 


(25) 


The  two  boundary  conditions  (15)  and  (l6)  serve  to  determine  the  coefficients 
a/TT^  and  of  the  two  exponential  terms,  each  of  which  is  a  solution  of  the  homo 

geneous  differential  equation  for  the  air  temperature.  The  following  two  equations 
may  be  solved  simultaneously  for  a  and  b: 


-TTjlel  -772 

e  a  +  e  b  |  + 


■] 


c  w 

a  I 

IhP  +C,+4C,T^ 

11 

r  -  -'^1 

uTTik-h  )e  a  +  (772k~h 

/)  b] 

-[ 

}p„Pc  +  (C,+4C,T^)  (P,+P,„)] 

k  [R  sin  ^le]  —  S  cos  <^|e|]  + 


hP  +  C|  +  4CjT 


k^t^R )  cos  co\  e\  + 


h  [hP„Pc  +  (Pc+Pm)] 

-(h'R+k^^S)  sin  a)|e|  tl-a]  hP^Q'  [^-h '  cos  <D|e|  +  uk  sin  cotefj 


(26) 


r  77, L  77, l1 

[_e  a  +  e  bj 


[■ 


77, L  77, L 

C_W  I  hP^  +  C,  +  4CjT"  11  77, e  a  +  77, e  b 


[hP„Pc  +  (C,+4C,T^)  (P^+P„)] 


[s  cos  coh  +  R  sin  col] 


cjC  W 
a 

+  C,  +  4C2TjJ^[^R  COS  COL  -  S 

sin  coL  J 

h 

’hP^Pc  +  (C,+4C,T^)  (P^+P„)_ 

+  [l-a]  hP^Q^Ct)  sin  coL 


Q'  is  determined  from  the  total  power: 


P  =  20' 


d  2 


TTz 

COS  - 

2L 


z 

i  .  i 


COS 


(27) 


(28) 
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where 


f  =  distance  from  flnx  maximum  to  end  of  fuel  element 
=  distance  from  center  of  pile  to  edge  of  loaded  lattices 
s  =  transverse  extrapolation  distance. 


It  is  found  that 


77-^P 


16  LN  [1+s/L']^  l^sin  7T/2  (1+s/L')^  [sin  (77-f/2L)  +  sin  (7r|e|/2L)] 


(29) 


The  power  transfer  coefficients,  one  for  the  air  in  the  gap,  one  for  the  air  in  the 
channels,  are  found  from  equation  (10).  The  constants  and  are  found  from 

experiments,  and  T^  may  be  determined  from  the  simple  calculation  discussed 
earlier.  The  mass  rate  of  air  flow  is  found  from  equation  (9). 


The  solutions  are  as  follows: 


a  77, z  b  77,2  S  R 

T  -  —  (e  -1)  +  -  (e  -1)  +  —  sin  + -  (l~cos  cJz ) 

77  ,  772  CjJ  0^ 


(30) 


C.w  hP_  +  C,  +  4C,T 


[■ 


77.  z 


77  ,  z 


a  +  e 


b  +  S  cos  coz  +  R  si 


in  CJzJ 


h[hP„P^  +  (C,+4C,T;^)  (P^+P„)] 


~  Q^hP^  (1”<^)  cos  Ct^z 


(31) 


o 

CD 

[C|  +  4C,T^] 

[• 

77  ,  2  TT^z 

a+e  b+Scos  coz  +  R 

sin 

h 

[h^P„Pc  +  (C,+4C,^^)  (Pe+P„)] 

+  (l“^) 


(32) 


Numerical  Results  and  Conclusions 


The  purpose  of  this  report  is  to  study  the  effects  of  various  parameters  on  the 
heat  distribution  in  the  Brookhaven  reactor.  The  results  of  this  study  are  summarized 
in  tabular  form  in  Figure  11,  and  in  graphical  form  in  Figures  3-10,  and  12. 
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Each  of  the  graphical  Figures  3-10  consists  of  six  graphs.  Three  of  these  graphs 
(solid  lines)  relate  to  the  quantity  being  studied  and  the  remaining  three  (dashed  lines) 
are  presented  for  comparison.  These  latter  three  refer  to  standard  values  of  the 
parameters  and  are  based  on  the  simple  method  of  calculation.  In  each  set  of  three 
graphs,  one  describes  the  fuel  element  temperature,  a  second  illustrates  the  graphite 
moderator  temperature,  and  the  third  relates  to  the  air  temperature.  The  air  temper¬ 
ature  as  it  enters  a  channel  from  the  gap  is  taken  to  be  zero  in  all  cases.  The  values 
of  the  parameters  taken  as  standard  will  be  found  in  the  table  defining  our  symbols. 

For  those  parameters  altered  in  each  graphical  figure,  both  the  standard  value  and 
the  altered  value  are  quoted.  Unstated  values  are  standard. 

The  tabular  Figure  11  is  a  comparison  of  the  fuel  element,  moderator,  and  air 
temperatures  when  one  parameter  at  a  time  is  changed  to  their  temperatures  under 
standard  conditions  as  computed  from  the  more  elaborate  theory.  The  standard  and 
altered  values  of  the  parameter  changed  are  indicated  in  each  case.  The  temperatures 
of  the  fuel  elements  are  compared  at  three  points:  at  the  gap,  at  the  plenum  end  of  the 
pile,  and  at  the  maximum  (which  may  be  shifted  under  the  altered  conditions).  The 
graphite  temperatures  and  the  air  temperatures  are  compared  at  the  gap,  at  the 
middle  of  the  reactor,  and  at  the  plenum  end.  A  positive  sign  means  the  temperature 
under  altered  conditions  is  higher  than  under  standard  conditions. 

The  results  allow  a  few  generalizations: 

1.  The  air  temperature  can  be  computed  as  well  from  the  simple  theory  as  from 
the  elaborate  theory. 

2.  The  metal,  graphite,  and  air  are  equally  hot  at  the  plenum  end  of  a  channel. 

3.  The  graphite  and  air  have  no  maximum  temperature,  at  least  of  any  con¬ 
sequence. 

4.  For  about  30%  of  a  channel  length  from  the  plenum  end,  the  graphite  tempera¬ 
ture  is  essentially  constant. 

5.  The  maximum  metal  temperature  occurs  about  half-way  along  a  channel. 

6.  The  heat  transfer  coefficient  to  the  air  from  the  graphite  or  metal  and  the 
specific  heat  of  the  air  are  the  two  parameters  of  major  importance. 

7.  The  flux  distribution,  the  fraction  of  heat  generated  in  the  carbon  and  the  heat 
conducted  or  radiated  to  the  graphite  are  of  secondary  importance. 

8.  The  thermal  conductivity  of  graphite  and  the  heat  transfer  coefficient  to  the 
gap  are  relatively  unimportant,  except  for  the  graphite  temperature  at  the  gap  face. 

Finally,  Figure  12  is  a  comparison  of  theory  and  experiment  as  regards  fuel 
element  temperature  for  two  channels.  The  agreement  is  poor,  even  though  the 
parameters  are  as  representative  as  possible.  Several  reasons  for  the  poor  agree¬ 
ment  may  be  adduced: 

1.  We  assume  the  heat  transfer  coefficient  to  the  air  in  a  channel  to  be  constant. 
Actually,  the  turbulence  is  greater  near  the  gap,  so  the  heat  transfer  coefficient  is 
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larger,  and  the  metal  is  cooler  than  our  calculations  indicate.  This  explanation  is 
one  of  the  few  which  will  account  for  the  rapid  increase  in  fuel  element  temperature 
near  the  gap. 

2.  Errors  occur  in  the  measurement  of  temperatures.  Thermocouples  were 
mounted  on  the  aluminum  fins  of  the  cartridges  and  were  located  directly  in  the  air 
stream.  These  thermocouples  cause  local  perturbations  in  the  air  flow,  hence  in  the 
heat  transfer  coefficient  and  in  the  metal  temperature, 

3.  Errors  take,  place  in  the  measurement  of  power. 

4.  We  assumed  all  air  going  through  the  reactor  goes  through  the  channels j  10% 
may  leak  around.  The  total  air  flow  was  measured.  The  pile  was  orificed  for  692 
channels  in  a  double  cosine  manner.  Of  these,  112  of  the  peripheral  channels  were 
plugged  and  assumed  to  transmit  no  air.  Of  the  remaining  580  channels,  2.5  channels 
had  no  fuel  element  in  them  and  no  orifice.  We  assumed  advisedly  that  these  trans¬ 
mitted  twice  as  much  air  as  an  unorificed,  unplugged  channel  loaded  with  a  fuel 
cartridge.  Not  only  errors  in  the  air  flow  per  channel  occur  in  the  foregoing  con¬ 
siderations,  but  air  may  leak  around  the  pile  at  points  other  than  the  channels.  Evi¬ 
dence  for  this  leakage  is  provided  by  the  unexpectedly  low  graphite  temperature.  Any 
leakage  would  affect  the  two  important  parameters  of  air  flow  per  channel  and  the  heat 
transfer  coefficient.  A  10%  decrease  in  air  flow  would  affect  the  maximum  metal 
temperature  by  roughly  20^0. 

In  the  calculations  relating  to  the  comparison  with  experiment,  we  assumed  the 
reflector  savings  to  be  61  cm  on  each  side  of  the  reactor. 4  The  maximum  of  the  flux 
was  calculated^  to  be  80  cm  from  the  gap  center. 

In  conclusion,  we  would  like  to  point  out  a  few  simplifications: 

1.  If  the  numerical  coefficients  are  examined,  especially  the  exponential  factors 
in  equation  (26),  to  a  very  high  degree  of  approximation  we  may  put  the  coefficient 
of  a  equal  to  zero  and  solve  directly  for  b, 

2.  Likewise,  in  equation  (27)  we  may  take  the  coefficient  of  b  equal  to  zero  and 
solve  directly  for  a. 

3.  In  view  of  the  unimportance  of  the  heat  transfer  coefficient  h'to  the  gap  air, 
h'  may  be  set  equal  to  0,  unless  the  finer  details  of  the  graphite  temperature  are  de¬ 
sired  right  at  the  gap. 

4.  Because  the  thermal  conductivity  k  of  the  graphite  is  of  negligible  impor¬ 
tance,  equation  (26)  becomes  finally 
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1  + 


C^W77, 


[hP„  P  C,  t  4C,t’] 


[ 


('■Pc*''''.) 


1 


-  ! 


CaW^^ 

4=  -Ip  cirk  r,)\  \  —  C  t,  . 

[hP^  +  C,  +  4C,T^] 

R  cos  co\e\  -1-  S 

>  sin  Cl)  1  e|  J 

(hPc+hP„)] 

+  (  1-Ot,)  sin  <^|e| 


(33) 


5.  Likewise,  equation  (27)  becomes 


-77, L 


1  + 


f +  (C,+4C,T^)  (hP^+hP„)] 


-  I 


— R  sin  coh  —  S  cos  <^L 


[> 


][■ 


Ca"'*'  U'^m  +  C,  +  4CjT;J  [_F  cos  + 
[h^P„Pe  +  (C,+4C,-?^)  (hP^+hP^)] 


-  S  sin  coh  \  +  (l-CL)  OJQ 


hP^  sin 
m 


coh 


(34) 


6.  The  air  temperature  is  found  from  equation  (30).  We  note  the  first  term 

is  of  consequence  only  at  the  plenum  end  of  the  channel;  the  second  term  is  important 
only  at  the  gap  end  of  a  channel.  The  sine  term  is  by  far  the  most  significant,  being 
at  all  times  one  and  a  half  orders  of  magnitude  more  than  all  others  combined. 


7.  The  carbon  temperature,  T^,  is  found  from  equation  (31).  The  term  in  a 
is  important  only  at  the  plenum  end  of  a  channel;  that  in  b  is  important  only  at  the 
gap  end;  the  cosine  terms  are  considerably  more  important  than  the  others,  usually 
by  as  much  as  one  order  of  magnitude.  The  term  T^  and  the  remaining  expression 
for  T^  are  of  comparable  size. 

8.  The  metal  temperature,  T^^,  is  found  from  equation  (32).  As  before,  the  a 
term  is  important  only  at  the  plenum  end,  the  b  term  at  the  gap  end,  and  the  cosine 
term  dominates  all  others  usually  by  at  least  an  prder  of  magnitude.  The  term  in  T^^ 
and  the  remaining  expression  for  T^  are  of  comparable  magnitude. 
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Figure  5.  Effect  of  specific  heat  of  air.  BNL  Log  No  D-1836. 


rigure  6.  Effect  of  actual  flux  distribution.  BNL  Log  No.  D~1838. 


Figure  7,  Effect  of  heat  generated  in  carbon.  BNL  Log  No.  D-1839‘ 
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T„  SIMPLE- STANDARD  CASE 


Conditions:  Cj  =  .019  cal/oC-sec-cm  (Cj  =  .0095  cal/oC-sec -cm  is  standard) 
cal/oC-sec-cm  (4C2  =  .0086  cal/oc-sec-cm) 
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^ COMPLICATED -THIS  CASE 
^ - Tq  SIMPLE -STANDARD  CASE 
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Figure  8.  Effect  of  heat  radiated  or  conducted  to  carbon.  BNL  Log  No.  D-1840 


Conditions:  k  =  0.40  cal/sec-®C-cm  (0.17  cal/sec-°C-cm  is  standard) 
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Unit 

X 

40.6 

101.5 

cm 

y 

61,0 

142.2 

cm 

w 

1640 
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Ibs/hr 

h 

.0053 

.0049 

cal/cm^- 

sec-OC 

4C2 

.0051 

.0033 

cal/cm- 

sec-°C 

Conditions; 

a  =  10%  P=  1 6.8  megawatts 

N  =  580  channels,  Lj'  =  274.5  cm 
but  2.5  not  loaded  -  218.5  cm 
L  =  292  cm 
le|  =  76  cm 
f  =  259  cm 

All  other  conditions  standard. 
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Figure  12.  Comparison  of  experimental  metal  temperatures  with  complicated  theory, 

BNL  Log  No.  D-1844. 
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